Ab initio methods for finite temperature two-dimensional Bose gases 
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The stochastic Gross-Pitaevskii equation and modified Popov theory are shown to provide an ab 
initio description of finite temperature, weakly-interacting two-dimensional Bose gas experiments. 
Using modified Popov theory, a systematic approach is developed in which the momentum cut-off 
inherent to classical field methods is removed as a free parameter. This is shown to yield excellent 
agreement with the recent experiment of Hung et al. [Nature, 470 236 (2011)], verifying that the 
stochastic Gross-Pitaevskii equation captures the observed universality and scale-invariance. 



Introduction — Ultracold atomic gases represent ver- 
satile tools with which to investigate many-body quan- 
tum physics [l[ . A key tunable property of these systems 
is their effective dimensionality; increasing the trapping 
potential, used in experiments to confine and cool atoms, 
in one or two directions produces gases that are effec- 
tively two-dimensional (2D) or one-dimensional (ID), re- 
spectively. Low-dimensional geometries in turn lead to 
richer physics due to the enhanced importance of fluctu- 
ations which restrict the onset of long range order d Q . 

A number of experiments have been performed recently 
in order to examine the thermodynamic properties of 
weakly interacting 2D Bose gases |4|-[l3| . Due to the high 
precision now routinely attained in such experiments, a 
particularly powerful feature is their usefulness in accu- 
rately testing microscopic theories. In this respect, the 
2D Bose gas is interesting as fluctuations are typically im- 
portant over a broad critical region [LH [l5j , meaning the 
standard mean field approach to weakly interacting Bose 
gases is not well suited; this broadness was however ex- 
ploited by Hung et al. [13J to obtain a clear experimental 
observation of critical phenomena in ultracold atoms near 
the Berezinskii-Kosterlitz-Thouless (BKT) phase transi- 
tion HE OS. 

At equilibrium, Monte Carlo (MC) calculations have 
been successfully applied to the uniform 2D Bose gas 
El El El , notably yielding a microscopic prediction for 
the BKT transition point [19( . The harmonically trapped 
case has been studied using quantum MC [2fJ, various 
mean field theories [2ll - |25j and classical field calculations 
[lEHU. Classical field methods have the advantage of 
providing a time-dependent description of the gas, which 
makes them additionally applicable to systems away from 
equilibrium; nonetheless, a simple yet accurate mean- 
field theory is also highly desirable to avoid the need 
for more complicated methods in calculating equilibrium 
properties. In this work we perform a quantitative com- 
parison between two such complementary methods: (i) 
the modified Popov (MP) theory of Stoof and co-workers 
[29l l30l ]. which is straightforward to solve and provides 
quick access to equilibrium properties in all dimensions, 
and (ii) the stochastic Gross-Pitaevskii equation (SGPE) 
[3TI - l33| ] (see also [34 |) that has already successfully de- 
scribed several ID Bose gas experimen ts f35l. 36| and is 
also applicable in dynamical situations [32|, |37|- The ex- 



perimental in situ measurements of Hung et al. [1 31 ] offer 
an ideal test-bed for the theories we wish to consider, free 
of the complications associated with modelling expansion 
imaging, thus a primary motivation for this work is to 
demonstrate the applicability of the SGPE in capturing 
the universality and scale-invariance of 2D experiments. 

A further related aim is to systematically circumvent a 
standard limitation of classical field calculations, linked 
to the fact that the Bose gas represented in this way 
must obey classical statistics. In practical terms, this is- 
sue manifests as a sensitivity to the numerical grid (or, 
equivalent ly, momentum/energy cutoff) used in solving 
such models [HH^], which, as a free parameter, can un- 
dermine their use for ab initio studies. While calculations 
based upon classical lattice models have been linked back 
to their quantum counterparts in the context of Monte 
Carlo studies [H H(| , the situation is less clear for meth- 
ods based upon dynamical equilibration of classical fields: 
approaches to date include selecting a cutoff based on the 
assumption that the classical field temperature should 
match that of an ideal Bose gas with the same conden- 
sate fraction (4l| , specifying a minimum acceptable mode 
occupation, typically within an energy cutoff somewhere 
in the range gzr>n < E cut < ksT [3^ ]. and a high temper- 
ature semi-classical field method valid for fc^T > fi [42l ]. 
Motivated by the desire for a formulaic approach to cut- 
off choice that is valid also at temperatures ksT < //, and 
which also includes the effect of interactions, in this work 
we have devised and tested a strategy by which to specify 
an 'optimum' SGPE grid choice, based upon comparing 
the classical and quantum limits of MP theory. 

Methodology — The 2D stochastic Gross-Pitaevskii 
equation describes the Bose gas via a noisy field tjj(x.,i), 
which satisfies the equation of motion 
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where V(x) = muj 2 (x 2 +y 2 )/2 = rr 
potential in the more weakly confined x-y plane, (72D = 
V / 87r(a// 2 )S 2 /to = gh 2 /m is the 2D interaction strength 
(with a the s-wave scattering length) , and 77 is a complex 
Gaussian noise term, with correlations given by the rela- 
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Figure 1. Optimum grid choice: (a) Variation in g' 2 '(0) with 
energy cutoff for MPrj at two different temperatures. Hor- 
izontal lines indicate the MPbe results; the horizontal coor- 
dinate of the point at which the MPrj data intersects the 
relevant MPbe result identifies the optimum energy cutoff. 
Dashed lines indicate the results when above cutoff atoms 
are neglected, (b) (x) and (c) density profiles from MPbe 
(solid black) and MPrj for several cutoffs: /3E = 0.20 (dotted 
blue), w 0.77 (dashed brown) and = 1.50 (dot-dashed red). 



tion (Tj*(x,t)r](x',t')) = 2tvy(x,t)k B TS{x-x')S(t-t'). In 
Eq. (JT|) n a bovc denotes the density of atoms with momenta 
greater than the momentum cutoff due to the numerical 
grid used to solve the SGPE H|. 

Classical field methods arise following the observation 
that the Bose field operator may be accurately replaced 
in the Heisenberg equation of motion by a complex val- 
ued field, under the condition that system modes repre- 
sented in this way are highly occupied. Although this is 
an excellent approximation in many circumstances, the 
system is then found to obey Rayleigh- Jeans statistics 
at equilibrium. In arriving at Eq.([T]), a similar notion is 
embodied in moving to a classical fluctuation-dissipation 
theorem (see Eqs.(38)-(40) of Ref. H|). 

The numerical solution to Eq.([T]) introduces an ultra- 
violet momentum cutoff implemented here by the dis- 
cretization scheme chosen; this implies that the equilib- 
rium thermal state that is achieved for a given set of 
physical parameters can differ quite dramatically through 
variation of the grid spacing alone. There is then clearly 
some ambiguity as to which cutoff choice provides opti- 
mum agreement with nature. To address this point for 
2D Bose gases, we choose to make use of the MP theory. 

The MP theory [2^, [3(| was formulated specifically to 
describe low-dimensional Bose gases, taking into account 
the effects of phase fluctuations to all orders (29|, and 
has been found to agree well with the SGPE in previous 
studies [M El El- In Ref.[H, MP was compared to 
both Hartree-Fock (HF) and MC results for the 2D Bose 
gas. The MP predictions for the critical chemical poten- 



tial and density for the BKT transition were in excellent 
agreement with those obtained from MC simulations [l!| . 
In fact, the MP theory was found only to break down in 
the region where its equation of state reduces to that of 
HF, which occurs close to the BKT transition due to the 
mean field nature of this approach. A renormalization 
group analysis was subsequently shown to avoid this dis- 
continuity and to match smoothly to the MP equation of 
state for chemical potentials /i > 0.18k B T (23j . 

The central ob ject within the MP approach is the 
quasi-condensate [461 ] which, in the Thomas-Fermi ap- 
proximation, obeys <?2D(^qc + 2nt) = /Lt e g where the effec- 
tive chemical potential may be written as /ip g = fx— V(x). 
The thermal density is calculated from [23, [3(J| 

I ^ x f e k 1 g 2 un qc \ 

(2) 

with fiwk = -\/ e k(ek + 2.g2D"qc) and £k = h 2 k 2 /2m. 

Optimum cutoff selction: Classical vs. quantum statis- 
tics — In the usual formulation of MP, = N^ E = 
l/(exp(/3fiiJk) — 1), since the particles obey Bose-Einstein 
statistics, and the upper index in the sum in Eq. 
k max , may be straightforwardly taken to infinity. How- 
ever, setting TVk = N^ J = l/{(3hu^) instead leads to 
a MP result based upon Rayleigh- Jeans statistics which 
is therefore analogous to using the classical fluctuation- 
dissipation result in the SGPE (see e.g. Eqs. (38)- (40) of 
Ref. HI). We will use MP RJ and MP BE to denote the 
Raylcigh-Jeans and Bose-Einstein cases, respectively. 

Since the MPrj approach gives very good agreement 
with the SGPE for a given cutoff [47j , this allows us to ex- 
tract an optimum value for the SGPE grid spacing, as we 
can examine the accuracy of the classical approximation 
as k max is varied by directly comparing the MPrj and 
MPbe results. To measure this difference, we compare 
predictions for the normalized second-order correlation 
function, chosen because of the important role played by 
density fluctuations in quasicondensation. In MP theory, 
this is given by # (2) = (n 2 c + 4n qc n t + 2n 2 )/(n qc + n t ) 2 . 
The procedure we adopt then is to first calculate g^ 2 \x = 
0) using the Bose-Einstein distribution in Eq. ©. This 
value then serves as a target result which we aim to 
match using MPrj with all other parameters held con- 
stant. Once an optimum cutoff is selected in this way, 
it is then used to calculate the optimum numerical grid 
spacing for use in the SGPE simulations. 

This approach is illustrated in Fig. [TJ symbols indi- 
cating the g( 2 >(0) values obtained from MPbe for two 
temperatures as the energy cutoff, shown on the horizon- 
tal axis, is varied. The diagonal dashed lines show the 
results if we neglect atoms above the cutoff momentum 
in the MPrj calculations. The symbols asymptote to- 
wards these lines for energy cutoffs > k B T, as there are 
then very few atoms in modes above the cutoff energy. 
For each temperature, the optimum energy cutoff may be 
read off from Fig. [1] as the horizontal coordinate of the 
point at which the classical g^ 2 '(0) data (squares/circles) 
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Figure 2. Density (top) and density fluctuations (bottom) from the SGPE, modified Popov (Bose-Einstein) and Hartree-Fock 
calculations versus the experimental data of Ref. JT^]. The rightmost plots show the SGPE results capture the experimentally 
observed scale invariance found in Ref. [l3|], including across the transition region indicated by the grey shaded area. 



intersects the horizontal line indicating the Bose-Einstein 
result. The lower plots show g^ 2 \x) and density profiles 
for the Bose-Einstein calculation (solid black) and the 
Rayleigh- Jeans for several cutoffs: (3E = 0.20 (dotted 
blue), 0.77 (dashed brown) and 1.50 (dot-dashed red). 

For the ideal gas, an expression was derived in Ref. [H| 
for a cutoff leading to optimum agreement between quan- 
tum and classical statistics; in 2D, it was found that this 
occured at an energy (3E w 1.6. Using the method out- 
lined above we typically obtain a value around f)E ~ 
0.77, which varies slowly with temperature in the range 
considered, yet is clearly much lower than the ideal gas 
prediction, indicating the important role of interactions. 

Comparison to Experimental results — In order to 
validate our scheme, we now directly compare the results 
to in situ experimental data of Hung et al. [H[ . We begin 
with the MP b e method and fix the trap parameters and 
temperature to those quoted in [13| . We then vary the 
chemical potential until the required density is reached 
at the trap centre. This provides us with all relevant 
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Figure 3. Universal behaviour of shifted density profiles (n — 
n c )A^ B versus rescaled chemical potential /3(/x — fi c — V(r))/g. 
Indicative errorbars are shown for several SGPE data points 
which originate from the fit used to obtain /i c . 



physical parameters required to solve the SGPE, barring 
an optimum choice of grid spacing, which is extracted as 
described; since there are no free parameters, both the 
MP theory and SGPE provide an ab initio description of 
the experiment. 

Figure [2] shows the SGPE (with an optimum momen- 
tum cutoff) and MP results together with the experi- 
mental density profiles of Ref. |1 31 ] . Corresponding data 
for density fluctuations, obtained from the normalized 
second order correlation function g^ 2 \r) via (Sn' 2 (r)) = 
(r?( 2 )(r) _ l)(n(r)) 2 , are shown on the bottom row of 
Fig|2] By comparing also to HF J49|, the trend we ob- 
serve is that despite matching the experimental density 
profiles well, the HF results consistently predict density 
fluctuations that are too large, relative to those observed 
experimentally. This may be understood since energy 
reducing correlations which lead to the onset of quasi- 
condensation are not incorporated in this theory, in ac- 
cordance with the findings for highly elongated gases pre- 
sented in Ref. [50|. In contrast, the SGPE and MP data 
agree very well with the experimental findings, thereby 
confirming the utility of each in describing the physics 
of 2D Bose gases at finite temperatures. Moreover, these 
findings highlight that choosing an optimum grid spacing 
based upon a higher order field correlation such as g^ is 
more robust, since the density is a less sensitive measure 
of the system. 

An important feature of Ref. (l3j was the experimen- 
tal demonstration of scale invariance in a 2D Bose gas. 
Plotting their data against a scaled chemical potential, 
/3(/i — V(r)), Hung et al. showed that densities measured 
at several temperatures collapsed to a single curve when 
appropriately scaled to the thermal de Broglie wave- 
length, AdB- This was found to be true in both the ther- 
mal and superfluid regimes; the fluctuation region may 
be identified with the crossover from a thermal gas to 
a superfluid, where the system moves from an enhance- 
ment to a suppression of density fluctuations respectively, 
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Figure 4. Scaled density fluctuation, (5-h 2 ) = X^ B (Sn 2 ), ver- 
sus the scaled compressibility, k = X 2 k//3 from SGPE (black 
circles) and experiment (red diamonds). Dashed and solid 
lines show (Sn 2 ) = k and (8h 2 ) = k/3 respectively. 



as may be seen in the lower rightmost plot of Figure [2] 
An important difference between the MP and SGPE ap- 
proaches, therefore, is the discontinuity in the density 
predicted by the former. In contrast, a key strength of 
the SGPE lies in its applicability across this fluctuation 
region, as may be seen in the rightmost plots of Figure [21 
which show that it also demonstrates scale invariance in 
good agreement with the experimental data. 

To estimate the location of the BKT transition from 
the SGPE simulations, we employ the fitting function 
used by Hung et ai, which yields a critical chemical 
potential, and a corresponding critical density, n c . 
These parameters can then be used to uncover the univer- 



sal physics of 2D Bose gases, since close to the transition 
point, the shifted density (n — nAAjg should be a univer- 
sal function of (/a— fj, c ) /3 / g alone [l3U14| . The universality 
of 2D Bose gases may then be tested by comparing scaled 
densities from gases with different interaction strengths, 
as in Ref.[l3|]. As shown in Figure EH by doing so we 
find that the SGPE reproduces the expected universal 
behaviour demonstrated experimentally. Moreover, we 
compare to the classical field MC results of Ref. 14j . 
again finding good agreement. 

Finally, to examine more closely how well the exper- 
imentally observed correlation effects are captured by 
the SGPE, we compare the density fluctuations to the 
compressibility. The compressibility may be calculated 
from the experimental data and SGPE using the expres- 
sion k = d(n) I 'cfyioff, the results of which are shown in 
Figure EQ which once again closely follow the results of 
Ref. 

Conclusions - We have found excellent agreement 
between ab initio calculations based on the stochastic 
Gross-Pitacvskii and modified Popov theories with in situ 
experimental data from the two-dimensional Bose gas ex- 
periments of Hung et ai. Combining these two methods, 
we devised and tested a systematic approach to circum- 
vent the problem of cutoff choice within classical field 
simulations that could prove crucial in modelling out-of- 
cquilibrium scenarios, which forms a natural extension to 
this work. 
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